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Abstract 

qq ■ This paper gives a new, simple and concise derivation of brane actions and brane 

dynamics in general relativity and in Einstein-Gauss-Bonnet gravity. We present 
a unified treatment, applicable to timelike surface layers and spacelike transition 
■ layers, and including consideration of the more difficult lightlike case. 

in ' 
o 

1 Introduction 

a -1 
u 



X 



Thin walls and shells of matter have surfaced increasingly in a variety of situations in 
astrophysics, cosmology and quantum gravity. Highly compressed expanding shells of 
material emerge from supernova explosions [1 , as false- vacuum bubble walls in inflation- 
ary phase transitions j2j, and in hypothetical scenarios of "new-universe" creation [3]. In 
a sudden global phase transition, the transition region can sometimes be idealized as an 
infinitesimaly thin spacelike surface layer Theoretical exploration of basic issues of 
principle, such as the possible outcomes of a classical or quantum gravitational collapse 
are often simplified, for purposes of a first reconnaissance, by idealizing the collapsing 
matter as a thin shell, thus reducing the complex differential field equations to simple 
algebraic junction conditions. 

With the advent of brane- world scenarios |6j, the scope of thin-shell dynamics has 
broadened to embrace higher dimensions and string-inspired extensions of Einsteinian 
gravity, in particular Einstein- Gauss-Bonnet (EGB) theory jjj. Recently, a class of Weyl- 
conformally invariant p-brane theories, which includes lightlike branes, has been proposed 

IE!- 

For purely Einsteinian shells the classical dynamics is straightforward [Hj (though the 
variational and Hamiltonian aspects can be subtle JH])- Appending Gauss-Bonnet terms 
to the usual Einstein-Hilbert action, however, is attended by a considerable step-up in 
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complexity. Initially, there were even doubts whether EGB shells admit a distributional 
description at all: the EGB bulk field equations develop ill-defined products of delta- and 
step-functions in the thin-layer limit unless the terms are arranged with care [TT]. The 
key is to express the field equations in canonical form, with distance from the layer in 
the role of "time" ; this segregates the most singular terms into the canonical momentum. 
The canonical momentum is the dynamical variable that "jumps" cleanly at a thin layer. 

The key role of canonical momentum suggets that the EGB junction conditions are 
most easily derived from the action. This derivation was carried through by Davis [12 , 
so that the basic equations of EGB shell dynamics are now well-established and widely 
employed [13]. Still to be desired is a systematic, self-contained exposition which draws 
together general-relativistic and EGB shell dynamics within a unified Lagrangian frame- 
work, and includes consideration of the lightlike limit [T3|. We hope this paper will go 
some way toward filling this gap. 



2 Toy Model 

To illustrate the essential ideas, we take a simple example from one-dimensional particle 
mechanics. We choose a "bulk" action functional of the path q = q(t), 



Sbuikiq] = (L - V ext ) dt , (2.1) 



with an acceleration-dependent Lagrangian of the form 

L = -b{q)q- l -b\q)q 2 -V(q), (2.2) 

where the functions b, V are arbitrary and V ex t(q,t) is an arbitrary external potential. 

We are considering q as an analogue of the metric, t as an analogue of distance normal 
to a boundary surface or layer, and q as an analogue of extrinsic curvature K ~ dg/dn. 
The two terms of (j2.1j) simulate the geometrical and matter actions; the external force 
F = —dV ex t/dq is the analogue of material stress-energy T a p. The particular, quasilinear 
functional form ()2.2)1 is patterned after the Einstein-Hilbert (EH) and Gauss-Bonnet (GB) 
Lagrangians. (The EH Lagrangian is quasilinear in the narrow sense that the coefficients 
of the second-derivative terms are functions of the metric only, not its first derivatives. 
But the corresponding coefficients do depend on first derivatives in the case of the GB 
Lagrangian) . 

Although L involves second derivatives, its quasilinearity ensures that the Euler- 
Lagrange equation 

for the classical path is no higher than second order. The classical path is thus uniquely 
determined by fixing its two endpoints q^ qf. But this path does not extremize the bulk 
action ()2.1|) . because the endpoint velocities q^ qf can still be varied freely: we have 

SS bulk [q] = + f\ 8q(t) dt +[ P 5q- 6B(q, q) ]*f , (2.4) 
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where 

P = b'(q)q; B(q,q) = b(q) q . (2.5) 

To have an action that is extremized by the classical path, one must add a boundary 
term [TS] to the bulk action (EHJ) : 

S[q] = S Mk [q] + B( q ,q)\ll , (2.6) 

The extremal of the action ()2.6|) now depends solely on the endpoints g/ of the classical 
path, in accordance with the Hamilton- Jacobi equation 

5S extrem (qi,ti;q f ,t f ) = p 5q\\ f . , (2.7) 

More fundamentally, the boundary term B is needed to preserve the composition law 

S[l -> 2 -> 3] = 5[1 -> 2] + 5[2 -> 3] , (2.8) 

for an arbitrary continuous path joining the points 1,2,3 with a sharp bend at 2 (as used 
for example in the "zig-zag" definition of the path integral [16j). 
The Euler-Lagrange equation ()2.3|) can be expressed as 

dp dV . i , 

This is the analogue of the gravitational field equations, with the right-hand side rep- 
resenting the stress-energy of matter. To simulate a thin surface layer, we consider an 
impulsive force acting at time t$: 

F(q,t) = a(q)5(t-t ). (2.10) 

This will produce a discontinuity in the momentum p. Since all delta-function contribu- 
tions to the left-hand side are gathered into dp/dt, it is straightforward to integrate the 
equation of motion (|2.9|) to obtain the jump across the discontinuity: 

[\P\} = tr(q), (2.11) 

where = lim e ^ {p(^o + e ) — p(to — e )}- This is the analogue of the geometrical 

junction conditions at a surface layer. The key to deriving it is just the identification 
of the canonical momentum p - most easily from the Hamilton- Jacobi equation (|2.7|) . or 
from a suitable adaptation of Lagrange's definition applied to an equivalent first-order 
Lagrangian, see (|2.16|) below. 

The jump condition (j2.11|) would follow as 

dSi mp (q , gp) _ _ _ 

dq 

from an "impulsive action" Si mp , as a function of position q(t ) = q and pre- and post- 
shock velocities q(t ± e) = q^. This is given by 



S imp = - [B{q,q)]^+l - / V ext (q,t)dt, (2.13) 

J to— e 
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and it coincides with the bulk action (|2.1j) if we choose tf = to + e, ij = t — e. Alterna- 
tively, (J2.13)) is obtainable without integrating L through the shock, instead considering 
to — e and t + e as future and past endpoints respectively, each with its own boundary 
action B, so that the bulk action ()2.6|) becomes 

S[q] = + J^J Ldt V ext dt + + [B}H +e , (2.14) 

when an impulse acts at t . 

It should be noted that the close relationship ()2.13|) . (|2.14|) of the boundary action 
B to an impulsive action hinges on the special form (j2.2j) of the action, and does not 
extend to an arbitrary quasilinear Lagrangian. It does, nevertheless carry over to the EH 
action, and also to the GB action modulo removable K 2 dK/dn terms. The origins of 
this peculiar circumstances will emerge in Sec. 3. 

The bulk + boundary action (J2.6|) is really a thinly disguised first-order action: 

S[q] = f S {In - V ext ) dt; I x = \ b'(q) q 2 - V(q) , (2.15) 
JU 2 

with the standard definition 

dli 

' = ■«■ (2 - 16) 

Therewith everything relating to this mechanical model takes on a trivial appearance. 
Not so, however, for its gravitational counterparts: there, only the original, second-order 
Lagrangian is a geometrical object and a scalar; the split into a first-order Lagrangian and 
a pure divergence cannot be made in a co-ordinate invariant and boundary-independent 
way. One is essentially forced to retain the bulk + boundary formulation. 

Let us finally note that, because of the freedom to redefine the bulk part of the total 
action by adding a total derivative, the definitions of canonical momentum and boundary 
action are (trivially) arbitrary to the extent 

I 1 ^I 1 + j t f(q) ; B - B - f(q) ; p - p + f\q) . (2.17) 

This has no effect on the impulsive jump conditions 1)2. lip , because the arbitrary function 
f(q) is continuous at the shock. 



3 Brane Dynamics: Einstein-Hilbert Action 

We begin by focussing on general relativity and on higher- dimensional gravitational the- 
ories governed by the Einstein-Hilbert action. Our purpose is to derive the well-known 
junction conditions Pj, ^3] which determine the motion of a surface layer in such theories 
from the action. We shall present a unified treatment, applicable to timelike, spacelike or 
lightlike layers. 

In an (n + l)-dimensional spacetime, the Einstein-Hilbert bulk action is 

Sbuik = y k J £E H (g,r)d n+1 x, (3.1) 

where 

C EH = V=gR = yf^gsT Wt, - ) - A , (3.2) 
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Ci = J=ggT{T% r^-r^i*,) . (3.3) 

The second-order Lagrangian (|3.2|) is degenerate, and can be reduced to first-order form 
by extracting a pure divergence: 

C EH = ~d x a x + C x , (3.4) 



where 

which is reducible to 



a\ g ,r) = v^(^r; Q - <rr*„), (3.5) 
1 



a 



d,(-gg x n- (3-6) 



The first-order Lagrangian C\ dates back to Lorentz, Hilbert, Einstein, Weyl and Felix 
Klein, and was employed by Einstein to define his pseudo-tensor for gravitational 
energy and radiation. 

The complete Einstein-Hilbert action Seh complements the bulk action (j3.1|) with a 
term Sm v coming from the boundary, which soaks up the pure divergence in ()3.4j) [TT] . 
Including also the matter contribution, the complete action reads 



SeH — Sbulk + Sbdy + S mat — — 



l - J C l( i n+1 x + S mat . (3.7) 



Here 



S bdy [g,T) = i- j d x a x d n+l x, (3. 



is convertible to a surface integral over the boundary S, and S mat yields the material 
stress-energy tensor 

3£>mat 1 .... fl , 

-gT^. (3.9) 



8g^ 2 

Computing the variation of the bulk action (|3.1j) is facilitated by the Palatini identity 

5R, U = 5I* MX - ST^ V , (3.10) 

where | stands for the covariant derivative associated with the metric g^ v . Recalling (1331) . 
and momentarily treating g and T as independent, this gives 

5C EH (g,T) = (5 g + 5 T ){y^g^R^(r)} = Sg^y/^G^ - d x (5 r cx x ) . (3.11) 

From (|3.7|) - (|3.9|) we thus obtain 

SS EH = ttW {Gpis - KT tlu )6g^V=gd n+1 x + 5 9 S bdy [g,T}. (3.12) 

The last term depends only on the metric variation 6g a p at the boundary, not ST. Thus, 
for metric variations which vanish at the boundary, we find 

J^7 = ° G, v = kT, u , (3.13) 

which are the bulk equations. 
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While the second-order Lagrangian (|3.2|) is a geometrical invariant, its split (|3.4j) into 
a first-order part and a pure divergence is co-ordinate-dependent. This split can, however, 
be endowed with geometrical significance by a special choice of the bulk co-ordinates x^ 
which anchors them to the boundary E. In general, E will be characterized by parametric 
equations 

E : = x"(f a ,z) ; z(x f *)=0, (3.14) 
where £ a (a = 1, ...n) are arbitrary intrinsic co-ordinates. Then ()3.8|) integrates to 

S bdy = 7T / ^9 x zdA, (3.15) 

where dA is a (non- invariant) element of n-area. We now impose the anchoring condition 
on x^: 

E : x° = z = , x a = C- (3.16) 

Then (|3.15|) becomes 

Suy = ^ j^a Q {g,T)d\, (3.17) 

with 

a = a x d x z = {-g)-^d x [{~g)g^d,z\. (3.18) 

Its non-covariant appearance notwithstanding, (J3.18j) is actually (twice) the mean 
extrinsic curvature density of E, up to a dynamically irrelevant term. Moreover, it is 
regular even for a lightlike E, so that ()3.18|) may be considered a valid extension of the 
notion of mean extrinsic curvature density to the lightlike case. 

To verify these statements, suppose first that E is non-lightlike, and thus has a unit 
normal n M , which is transverse to E, say in the direction of increasing z, (n° > 0); n.n = e 
with e = +(— )1 whenever E is time(space)like. Then ()3.18|) can be shown to reduce to 

a = 2e/C + ^g~eN- 2 d a s a ; IC = \ {n) g\^ K . (3.19) 

where K = n a \ a is the mean extrinsic curvature, and N, s a the lapse and shift which 
appear in the standard ADM form of the bulk metric 

ds 2 = g ab (dC + s a dz) (d£ b + s b dz) + eN 2 dz 2 , (3.20) 

so that 

g 00 = eN -2. g 0a = _ eN -2 s a. |(n)^|| = ^ jy-1 _ (33^ 

In the anchored co-ordinates (J3.16|) . n M has components 

Uil = eNd^z, Nn° = e, N n a = - s a . (3.22) 
It is now straightforward to derive (ETPIJ) : from (l3~T£l) . (I3~2"T1) and (l3~221) . 

<x°-2e/C = {-g)--*d^{-g)eN~ l n»} - 2eN' 1 d fl (^n^) 

= -y/=geN- 2 d„(Nn») = y^g eN~ 2 d a s a . (3.23) 

This extra contribution ()3.23|) to the boundary action ()3.17|) , ()3.19|) is "velocity-independent" 
(i.e., independent of transverse derivatives dog^), and has the same role as /(g) in ()2.17|) . 
It could be reconverted to a volume integral and included with L\ in the first-order bulk 
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action in (J3.7|) . Its contribution to the momentum is trivial. Moreover, it contributes 
nothing at all to the jump of momentum across an extrinsic-curvature discontinuity, 
since its jump across E is zero if x M and g^ are continuous. This will continue to hold in 
the lightlike limit (N — > oo). 

Recalling (|3.5j) . we conclude that 

e | JC(g, r) | = C?(g) [| |] ; C?(g) = ^ (g^ - g^Sf) d a z , (3.24) 

represents precisely the jump of mean extrinsic curvature density at an arbitrary (lightlike 
or non-lightlike) extrinsic-curvature discontinuity S. (For lightlike E, e is defined by 
continuity with the timelike case.) 

A surface layer can be characterized by supposing the bulk divided into two subdo- 
mains z < and z > with edges £_(z = —0) and X+(z = +0), glued together to form 
a common boundary E, loaded with a surface distribution of stress-energy; S + and E_ 
are supposed intrinsically isometric, but d a g^ u undergoes a jump at £ in smooth (e.g. 
skew-Gaussian) co-ordinates, i.e., £ is an extrinsic-curvature discontinuity. Dynamically, 
the surface layer is accounted for by adding to the Einstein-Hilbert action (J3.7|) a shell 
contribution equal to the sum of the actions (|3.17j) for the boundaries £_ and S + of the 
subdomains, taking into account the opposite directions of their outward normals: 

SEH,shell = 7^— / 2e[|/C|] d n £ + S m at, shell ■ (3.25) 

The second (matter) term S matjS heii generates the surface stress-energy density S af3 via 



SS, 



mat, shell 



Sg, 



af3 



- S al3 , (3.26) 



in analogy with ()3.9|) . 

Variation of the total action (|3.7|) + (j3.25j) then yields the jump conditions 

[|^|] = \s a *>, (3.27) 

in addition to the bulk field equations ()3.13|) . 

The canonical field momentum density 7i at3 associated with any boundary S is defined 
in anchored co-ordinates (j3.14j) by 



p 1 dCx 

7T H ' 



(3.28) 



2k dg a p :0 

It is more conveniently extracted from the Hamilton- Jacobi variational formula 

8S EH = J ^ Sg a p d n i = 5 g S bdy (g, r) , (3.29) 

modulo the bulk field equations (|3.13|) . and recalling (J3.12)) . For the shell (|3.25|) . we obtain 
from (HEUD, (I3~T7|1 and (l3~Tn|) the jump 

K 7T p U = , (3.30) 

O9op 
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with g, T treated as independent, in analogy with the bulk identity 

V^G a , = ^_ J V^Rd n+1 x = (v^^^(r)) . (3.31) 

Explicit evaluation of ()3.30|) with the help of ()3.24j) yields 

[I ^° |]=0, (3.32) 



MK b |] = \v=9l\9c d ,o\]{g 00 (g ab g cd - g ac g bd ) 

+ 2s°V )( y )0 - g a0 g bo g cd - g ab g c0 g d0 }- (3.33) 

This is equivalent to a result (eq.(17) of Barrabes-Israel [13]) previously derived by inte- 
gration of the field equations through the layer. 

Like the Bianchi identity G al3 \/3 = for the bulk field equations (j3.13J) . the transver- 
sality condition (J3.32|) . 

[| *" p I] {d p z) = , (3.34) 

may be regarded as a consequence of the co-ordinate-invariance of the action. The bound- 
ary action f)3. 17j) is invariant under the infinitesimal anchored co-ordinate transformation 

x" -> = x^ + ^2?rf(x). (3.35) 

This gives 

hg a /3 = g a p(x) - g a p{x) = 2r] ia d p) z + 0{z) , (z -»• 0). (3.36) 
Hence from (|3.29p . 

6lSeh = J ^ 2 Va dpz d n i , (3.37) 

which is required to vanish by co-ordinate invariance of the total action, leading to (|3.33|) . 
For a non-lightlike layer ((? 00 |s 7^ 0), (|3.33|) simplifies to 

« I * ab \\ = \V=g[\ g cd ,o |] g°° (A ab A cd - A ac A bd ) (3.38) 

where A ab projects onto S and coincides with the inverse intrinsic metric in anchored 
co-ordinates 

„a0 „60 

A ab = g ab _ g_g_ = g ab _ ^arf = (n) g ab _ (3.39) 

The jump conditions ()3.27)1 can then be reduced to their standard non-lightlike form [Oj 

- ^\ {n) g\ 1 *[\K ab - {n) g ab K\] = k [| 7r ab |] = ^\ {n) g\?S ab , (3.40) 
in terms of the jump of extrinsic curvature 

n^o = ^D^n = ^i^d^d, (3-41) 

and of the surface stress-energy tensor S of the shell 
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4 Gauss-Bonnet Action 



When the Einstein-Hilbert action (|3.1|) is augmented with term quadratic in the curvature 
the simple form (|3.4U|) of the junction conditions is no longer valid. In fact, a distributional 
brane dynamics is no longer even possible in general, because the bulk field equations now 
involve inadmissible products of distributions in the thin-layer limit. The exception is the 
case when the quadratic terms have the Gauss-Bonnet form. In this case the bulk field 
equations are quasi-linear, and a distributional description of thin layers remains viable. 
In this Section, we shall examine how the junction conditions (|3.40j) are modified. 
The Gauss-Bonnet action is (see the Appendix for further details) 

Stuik = £ / C GB (g,T)d n+1 x, (4.1) 



2k 

where a is the Gauss-Bonnet coupling constant, and 

£gb(9, T) = V^gK = J=j ~ 9 25 9 2 ' 5 ' #684 R V vv« • (4-2) 

Because of the plethora of indices in such expressions, we have found it convenient 
in many instances to let numerical indices 1,2,..., l',2', .... do duty for literal indices 
«i, «2, «2, •••• They are understood to run from to n, and repeated indices are to 

be summed. In contrast, we reserve the index to stand just for its numerical self, and, 
as in (jHHEJ), x° = z = will represent the boundary E. 

There is no holonomic split, analogous to (|3.4p . of Cqb into a first-order piece and a 
pure divergence. We must therefore proceed more indirectly to find the supplementary 
boundary action which effectively removes second derivatives from the GB bulk action 

Varying the affine connection V in (J4.2j) . and noting the Palatini and Bianchi identities, 
5-R 1 534 = — 2 5r53| 4 and <5; 43 4 R 1 5/3/4/14 = 0, we see that the T-variation of Sbulk involves 
a pure divergence, convertible to a surface integral: 



— SrS M k = I s/=g~(T X \xd n+1 x = I ^ga A d x zd n x, (4.3) 



a 

where 



a 



S^g^R^^ST^. (4.4) 

This is in complete analogy with (|3.11|) in the Einstein-Hilbert case, and it means that 
the bulk field equations 

l_ d£(g,T) 1 SS mat = , , 

2k dg^ 2 5g^ ' 1 ' ' 

with C = Ceh + oiCgbi are obtainable by simply differentiating the bulk Lagrangians 
(14. 2|) and (|3.2|) partially with respect to g. 

When the bulk field equations (|4.5|) are satisfied, the boundary term (|4.3|) gives the 
total variation of the bulk action. This involves ST, from which the variations Sg a pfl 
of transverse derivatives -i.e., of extrinsic curvature in the non- lightlike case - must be 
removed by compensating variation of a suitable boundary action. 

To isolate these extrinsic curvature variations, we assume for simplicity that S is 
non-lightlike and introduce Gaussian co-ordinates based on E as in (|3.2(J|) . Then 

(n+1) r^ = {n) T c ah (n+1) C = eK ab , ^T a ob = K\, (4.6) 
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where the extrinsic curvature K ab = ~ d g ab , and Latin indices run from 1 to n. We note 
also the Gauss-Codazzi relations 

{n+1) Ra bcd = {n) R ab cd - 2eK a[c K d]b} (4.7) 

{n+1) R^c d n» = 2K b[c . d] , (4.8) 

where ; represents the covariant derivative associated with the n-dimensional metric g ab . 
Retaining only the variations 8K ab in (|4.4jl - i.e., assuming 5g a p\x = - we find 

a = -5^,R^'^2e6K\, (4.9) 

so that (|4.3jl becomes 

-SkSmh = -2e5™'4 I ~ 2eK v v K*«) 5K\y/=g<Tx, (4.10) 

where we have made use of (|4.7|) . 

The boundary action S bdy must be chosen so that its if- variation cancels ()4.10|) : 

5K(S bu ik + S bdy ) = 0, (4.11) 

Since the intrinsic Riemann tensor ^R is independent of K ab , the choice 

^S bdy = 2e6° ff4 / s (W' 2 '^ - j&'yK*^ K^d^x, (4.12) 

meets this requirement, a result originally due to Myers [T5] . 

We are now ready to derive the Gauss-Bonnet field momentum 7r ab associated with E 
from the Hamilton- Jacobi equation 

S(S bulk + S bdy ) = J n ab 5g ab dJ:, (4.13) 

modulo the field equations ()4.5p . 

Evaluation of the left-hand side of ()4.13|) is simplified by noting that there is no K- 
contribution because of 1)4. lip , and none from the metric factors in ()4.2p because of (|4.5p . 
So 5S bu ik = SrSbuik is again given by ()4.3j) and ()4.4j) . with oT now effectively determined 
solely by the metric variation Sg ab , with K ab fixed. We re-express ()4.4p as 

~0 _ 9 r303'4' 25 dI'O rpl _ a , s:03'34 k \ 2'5' Jfnl' f A 1 /h 

a — -zo 12VQ g K 3/4' 01 53 — 4e0 011 , 2 , A 3.4/^ 01 5/3/, l 4 - i4 J 

where we have interchanged primed and unprimed indices in the second line and made 
use of the Gauss-Codazzi relation ()4.8j) . 

Turning now to 5S b d y in f)4. 13j) . one piece of this arises from variation of the intrinsic 
curvature (n) i? in (|4.12j) : 

-5S bdy 2eo« I g^'5^R l '^K\\^g\U n x. (4.15) 

Applying the intrinsic Palatini identity S^R 1 5/3/4/ = — 2 oT^,.^, one sees that the 
integrands ()4.15|) and (J4.14)) add up to an intrinsic divergence, which may be discarded. 



10 



All that remains to account for is the metric variations arising from \J—g, the raised 
indices (2' on {n) R\., l'2'l on the ^-factors) in and from dT 1 ^ = 5g lb K b3 in the 

variation of the bulk action (4.3), (4.4). The result of a straightforward calculation is 

^rr ab = \^g\ 1 i{6K a[m K™KZ ] + 6K b[m K™K: } 

- Ag ah K{ l K™K n n] + AeK cd *^R\ cbd ). (4.16) 

where *^R* ac bd - see the equation (jA.7|) of the Appendix for its definition - is the left and 
right dual of the intrinsic curvature tensor of S. This is equivalent to results previously 
obtained by Davis and others [T2*] . 

Following the argument leading to ([3.27)1 . we conclude that the dynamics of a non- 
lightlike shell in Einstein-Gauss-Bonnet theory, with bulk action 

S b uik = ^ I V^(R + aU)d n+1 x + S mat , (4.17) 
is governed by the junction conditions 

D ** | = S ab ; 7r ab = ** B + ^ b B , (4-18) 

where S ab is the surface stress-energy tensor of the shell, and is defined by ()3.26|) and 
(|3.42|) . The momenta tt^ h and iXq B are given by the first of (|3.4U|) and (|4.16j) respectively, 
and jump together with the extrinsic curvature across the shell. The action due to the 
shell augments ()4.17|) with a surface term -cf ()3.25|) and (|4.12j) - and is equal to 

S sheU = ^ / (| B EH |] + a [| B GB \})dZ + S mat , shell , (4.19) 

where 

B EH = 2eK (4.20) 

Bgb = 2e5% ^ n) R ab cd - jK a c K b d ^J K) . (4.21) 

These results hold for non-light-like shells, for which extrinsic curvature is well-defined. 
We now add some remarks on the lightlike case. Since this is a special limit of the timelike 
case, one might at first sight expect the resulting junction conditions to be simpler. 
Actually, however, this is far from being the case. Lightlike discontinuities propagate 
along characteristics of the field equations. It is a nontrivial matter to disentangle the 
lightlike discontinuities due to the shell from the accompanying gravitational shock waves. 
For pure Einstein theory this is still quite manageable, and we have presented the results 
in Sec. 3. But we have not yet succeeded in reducing the lightlike junction conditions for 
Einstein-Gauss-Bonnet theory to a form that we consider worth publishing. The nature 
of the characteristics themselves is made more complicated by the fact that the field 
equations are now quasiliniear only in the broad sense (linear in second derivatives but 
with coefficients depending on first derivatives). 

It is, however, straightforward to obtain the Gauss-Bonnet boundary and shell actions 
in the lightlike case. Under an arbitrary metric variation, the total variation of the bulk 
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action, is still given (modulo the bulk field equations) by (|4.3|) and (|4.4j) . which we can 
express in the form 

a = -4*ir 3012 5i\ 2 3, (4.22) 

where we have introduced the left and right dual of the curvature tensor, defined by 
(jA.3|) and (jA.6|) in the Appendix. In (J4.22|) we must separate out (and neutralize with a 
boundary action) the contribution of variations 5g a /3 } o of transverse derivatives: 

8F[i,2]3 = -d[ 1 Sg 2 ] 3 -> - (d[iz)5g 2 ]3fl. (4.23) 
Hence the boundary action must have the compensating variation 

^ 5S bdy = 4 J *R* 0203 5g 23fi ^ cTx . (4.24) 

Now, when the surfaces, x° = const, are lightlike, the components *R* 0203 do not contain 
second transverse derivatives g a /3,oo and they are linear in first derivatives, i.e., 

*r>*0a0b TS&b , r^abcd „ jsabcd r^cdab /a nr\ 

R = K (o) + K (i) 9cd,o, = A(i) , (4.25) 

where Km, Km are independent of g ab $ {Km has a fairly complicated linear dependence, 
not reproduced here, on the "nominal extrinsic curvature" JC ab = — T ab , which depends 
only on the intrinsic geometry for a lightlike S ). From ()4.24j) and (|4.25J) we infer 

^S bdy = 4 J (k$ + l -K$f 9cd ^ g abfi V=9d n x, (4.26) 

as the form of the Gauss-Bonnet boundary action. 



5 Concluding Remarks 

We have given an elementary, self-contained derivation of the action (useful for calculation 
of quantum tunneling amplitudes) and dynamical equations (i.e. junction conditions) for 
thin shells and branes in Einstein-Gauss-Bonnet theory. Our exposition has attempted to 
integrate as far as possible the treatment of timelike, spacelike and lightlike layers. For 
lightlike shells, the dynamics is complicated (especially in the Gauss-Bonnet case) by the 
fact that gravitational shock waves will in general accompany the shell. However, this 
problem should be ameliorated in situations of high symmetry, and this is currently under 
investigation. 

A Notations for GB 

In 1932, Lanczos ^H] noted that the Lagrangian 

^ = ^ ^121'2' -R 12 34 -R 1 2 3'4' , (A.l) 

leads, like the Einstein-Hilbert Lagrangian 

R = ~<^ 12 34, (A.2) 
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to field equations which involve no higher than second derivatives of the matric. R and 
1Z are the first two members of a family of Lagrangians having the same property found 
by Lovelock [20] ■ The n th member involves a product of n curvature factors formed by 
an obvious generalization of (jA.l|) . (jA.2|) . It has the "Gauss-Bonnet" property of being a 
pure divergence in a space of dimension 2n, and it vanishes identically in spaces of lower 
dimension. Properties of the Lovelock family are reviewed by Meissner and Olechowski 
[7j and Deruelle and Madore 

By defining the left and right dual of the curvature tensor 

*p* 34 _ 1 r343'4' pl'2' / A o\ 

-K 12 — -g 0121'2' M 3'4' , l A -<jJ 

and using the identity 

5343'4' = 4!4 3 5 2 4 44' 1 

= 5? 2 5U + 4 4 ' - 2 ST 5[t' ] - 2 5f , (A.4) 

where 5f| = 8f 8 2 — 5f 5 2 , f)A. If) and (|A.3|) can be recast as 

U = R 12 U *R\ 2 U = (R 12U ) 2 - A(R 12 ) 2 + R 2 , (A.5) 

(which was the form originally given by Lanczos). The tensor *i?* 12 34 generalizes to 
spacetime with dimension higher than four the left and right dual of the curvature tensor 
of general relativity It is equal to 

*iT 12 34 = i? 12 34 - 4 5[ 3 R% + ~ 4 4 R = tf 12 34 - 4 5[ 3 G\ - \ 511 R . (A.6) 

A similar quantity *^R* a b cd can be defined for the intrinsic curvature tensor 

*(n) R * ab cd = }_ 6 cdc> /d > (n) R a'b' cid/ = (n) j^cd _ 4 g[c { n) R d] ^ + ^cd (n) R ^ (A J) 

where the latin indices run from 1 to n. 

Because, as noted in Sec. 4, the T-variations of R and 1Z are pure divergences, the 
bulk field equations 

Gfa, + aQfa, = kT^, (A. 8) 

can be obtained from the bulk action 



Sbulk — 2 



r [V=g(R + aTZ) + S mat , (A.9) 



simply by partially differentiating the Lagrangian, holding i? 1 234(F) fixed: 



1 d 1 
1^ = -=—(^/=g~1l) = 27^ - n9r>K, (A.10) 



where 



-g dg^ v v 3 ' ^2 



n = g^TZ^; = -g ap (TZ aflPl/ + TZ au ^) , (A. 11) 

R\v2<- (A.12) 
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(jA.10|) confirms that the field equations follow the action in containing no higher than 
second derivatives of the metric. Noether's theorem and co-ordinate-invariance of R and 
1Z imply the contracted Bianchi identities 

= Q aP \p = 0, (A.13) 

which ensure compatibility of (jA.8|) with the conservation law T al3 \p = 0. All the above 
considerations extend straightforwardly to higher members of the Lovelock family. 
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